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Dissipative dark matter, such as mirror dark matter and related hidden sector dark 
matter candidates, require an energy source to stabilize dark matter halos in spiral 
galaxies. It has been proposed previously that supernovae could be the source of this 
energy. Recently, it has been argued that this mechanism might explain two galactic 
t^. | scaling relations inferred from observations of spiral galaxies. One of which is that 

p r is roughly constant, and another relates the galactic luminosity to r . [po is the 
dark matter central density and r is the core radius.] Here we derive equations for 
the heating of the halo via supernova energy, and the cooling of the halo via thermal 
bremsstrahlung. These equations are numerically solved to obtain constraints on the 
Po, ?"o parameters appropriate for spiral galaxies. These constraints are in remarkable 
agreement with the aforementioned scaling relations. We examine also constraints 
on the dark matter halo parameters of dwarf spheroidal galaxies. As with spirals, 
the dark matter structure of dwarf spheroidal galaxies can be similarly explained. 
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1 Introduction 

The standard ACDM scenariofl], which invokes weakly interacting dark matter par- 
ticles, has proven to be extremely successful in explaining the observed large scale 
structure and cosmic microwave background anisotropics [2]. However, this scenario is 
challenged by observations on smaller scales. For example, measurements of rotation 
curves of spiral galaxies indicate the existence of a dark matter core profile (e.g. [3]) 
in disagreement with the cuspy profile predicted by simulations of non-interacting 
dark matter [4]. 

Another challenge to the standard ACDM scenario is the lack of any evidence for 
new stable particles in LHC collider searches[5]. From a particle physics perspective, 
though, there is no compelling reason to favour weakly interacting particles over 
richer dark matter scenarios. In fact, the success of the standard model in explaining 
all collider data to date is a definite hint that dark matter resides in a hidden sector. 
That is, the fundamental Lagrangian decomposes into two sectors, one describing 
the standard particles and forces, and another which will contain the dark matter: 

£ = £-SM + C-dark ■ (1) 

The sector describing the ordinary particles has G = SU(3) <S>SU(2)l<S>U(1)y gauge 
symmetry, while the dark sector has independent gauge symmetries, G' . In this case 
one expects dark matter to be multi-component, self- interacting and perhaps dissi- 
pative, if the hidden sector has an unbroken U(l)' gauge symmetry. Such scenarios 
can also successfully explain the large scale structure and CMB, but can yield very 
different physics on small scales. Such hidden sector dark matter has been discussed 
in the literature in a variety of contexts, see. e.g. [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 
17, 18, 19, 20]. 

If the hidden sector does indeed contain an unbroken U(l)' gauge symmetry, then 
the associated 'dark' photon can kinetically mix with the ordinary photon: 

L-mix — yf r ^ v [2) 

where F^ u (F' ) is the field strength tensor for the photon (dark photon). Such 
kinetic mixing is renormalizable[21] and e can be viewed as a fundamental parameter 
of the theory. The physical effect of the kinetic mixing interaction [22] is to induce a 
tiny ordinary electric charge (oc e) for the hidden sector U(l)' charged particles. 

Mirror dark matter is a well motivated and concrete example of hidden sector dark 
matter in which C^ark is an exact duplicate of the standard model sector [23]. This 
means that the hidden sector has gauge symmetry G' = S , £/(3)'(g)S'£/(2) / L cg)£/(l)y. If 
the chiral left and right handed fermion fields are interchanged in such a hidden sector 
then the theory has an unbroken 'mirror' symmetry mapping each ordinary particle 
onto a 'mirror' particle (along with x —¥ —x). That is, for each ordinary particle, 
e, is, u, d, ...,7, ... there is a corresponding mirror particle, which we denote with a 
prime ('): e',u',u',d', ...,'-/', .... The unbroken mirror symmetry ensures that the 
mirror particles have the same masses as their corresponding ordinary counterparts. 



Similarly the gauge self-interactions (mirror electromagnetism etc) have the same 
form and strength in the mirror sector as they do in the ordinary sector. Although 
our focus here is on mirror dark matter, our results may be relevant to closely related 
hidden sector models, such as the ones discussed in ref.[15]. 

Mirror dark matter has emerged as an interesting dark matter candidate, for re- 
views and more complete bibliography see e.g. [8]. Mirror dark matter can explain[9] 
the large scale structure of the Universe - the matter power spectrum and CMB - 
in a manner analogous to standard collisionless cold dark matter models provided 
that e ~ 10~ 9 [24]. Such values of e are well consistent with direct laboratory lim- 
its, which arise from rare decays of orthopositronium[25]. Importantly, mirror dark 
matter can also explain[26] the positive dark matter signals from the DAMA[27], 
CoGeNT[28], CRESST-II[29] and CDMS/Si[30] direct detection experiments. This 
requires photon-mirror photon kinetic mixing of strength e ~ 10~ 9 . Of course, this 
type of rich dark matter candidate can feature a whole range of phenomena, espe- 
cially on small scales (galactic and smaller) which have only begun to be explored 
(see e.g. [31, 32, 33, 34, 35]). 

In this article we return to the problem of galaxy structure in the context of 
mirror dark matter. The existence of the unbroken U(l)' interaction will inevitably 
lead to significant self-interactions of the mirror particles. An implication of this is 
that galactic halos of spiral galaxies would have to be composed (predominately) of 
mirror particles in a pressure supported spherical plasma[33]. There may also be 
a subcomponent consisting of compact objects such as old mirror stars. Some of 
the interactions within the halo will be dissipative, such as thermal bremsstrahlung 
(e.g. e' + He' — >■ e' + He' + 7') which can cool the halo. At first sight, this might 
put into question the very existence of the halo. At any rate, the stability of the 
halo needs to be explained. A possible explanation was suggested sometime ago [33]. 
The idea is that ordinary core collapse supernovae provide the required energy. In 
the hot and dense core of a type II supernova mirror electrons and positrons can 
be created from kinetic mixing induced plasmon decay processes [36]. Thus ordinary 
supernovae can be a source of light mirror particles as well as the ordinary neutrinos. 
Indeed, it is estimated that mirror particles (e', e', 7') carry off roughly half of the 
core collapse supernova energy if e ~ 10~ 9 [36, 34]. A significant fraction of this 
energy can be absorbed by the mirror particle halo and thereby potentially replace 
the energy lost due to dissipative interactions. Order of magnitude estimates [33] 
suggested that the amount of heat generated roughly matched the energy dissipated 
for the Milk Way. The same mechanism could be responsible for stabilising the 
halos in all spiral galaxies. If so, then matching the heat supplied from ordinary 
supernovae to the energy dissipated suggests [37] a rough galactic scaling relation for 
spirals: Rsn oc Po r o- Here, po,r are the dark matter central mass density and core 
radius and R SN is the galactic supernova rate. It was further suggested that the 
energy from supernovae can predominately heat the inner regions of galaxies and 
that this would lead to a cored dark matter distribution. A second scaling relation, 
p r ps constant might be explained in this way, it was argued[37]. 

The purpose of this article is to provide a more detailed numerical analysis of this 



whole picture. Let us mention at the outset that the problem of small scale structure 
is a complicated one and the analysis performed here, although progress over the 
rough physical arguments of ref.[33, 37], still has some deficiencies. A number of 
assumptions are made, which would require further checks and refinements. More 
importantly, we do not attempt to evolve the galaxy from an early time to its present 
state, but rather see if we can at least explain some of the current properties of spiral 
galaxies. 

This paper is structured as follows. In section 2 we give a brief overview of the 
central idea, that dissipative dark matter candidates like mirror dark matter (and 
by extension, closely related hidden sector models), can have their halo's stabilized 
via supernova energy. In section 3 we discuss the hydrostatic equilibrium condition, 
which we numerically solve to obtain the galactic temperature profile for several 
example spiral galaxies. In section 4 we consider the ionization state of the halo. The 
equations governing the ionization fractions of the mirror helium, mirror hydrogen 
and mirror metal components are given and numerically solved. In section 5 we 
derive equations for the heating of the halo via supernova energy, and the cooling 
of the halo via thermal bremsstrahlung. These equations are numerically solved to 
obtain constraints on the po? r o parameters appropriate for spiral galaxies. These 
constraints are then compared with the scaling relations inferred from observations 
of galactic rotation curves. In section 6, we apply the same ideas to dwarf spheroidal 
galaxies and in section 7 we make a few comments on elliptical galaxies and galaxy 
clusters. Finally we give a few concluding remarks in section 8. 

2 The heating of the galactic halo 

The physical picture is that spiral galaxies such as the Milky Way are currently 
composed of ordinary matter in a disk, and mirror dark matter predominately in 
a (roughly) spherical halo. The halo consists of an ionized plasma formed out of 
the mirror particles, e',H',He',0',Fe', .... The plasma dissipates energy due to 
thermal bremsstrahlung and other processes and this energy needs to be replaced. 
The idea [33] is that ordinary supernovae can supply this required energy if photon- 
mirror photon kinetic mixing exists, Eq.(2). Such kinetic mixing gives the mirror 
electron and positron a tiny ordinary electric charge of magnitude ee. The energy 
loss rate for production of such minicharged particles from supernovae has been 
estimated in Ref. [36]: 



Qp = ^4e V [£ + ^- T*Q X (3) 




where Q± is a factor of order unity, and \x e is the electron chemical potential and 
T sa 30 MeV is the temperature of the supernova core. Demanding that Qp does 
not exceed the energy loss rate due to neutrino emission implies that e ~ 10~ 9 [36]. 
Thus, supernova can be a source of energetic light mirror particles e', e', 7' which 
can ultimately replace the energy lost due to radiative cooling. This heating of the 



halo is in the central region of the galaxy, which leads to the temperature having a 
mild radial dependence. The halo is generally hotter at the center and decreases as 
the distance from the center, r, increases. The 'average' halo temperature (say at a 
distance r = 3ro where r^ is the disk scale length) is typically of order 300 eV for 
the Milk Way, and ranges from 10 eV for the smallest spirals to around few keV for 
the largest spirals (see the discussion in the following section). 

The amount of supernova energy required to replace the halo energy lost due to 
radiative cooling is sizable. Estimates[33, 37] indicate that at least a few percent 
of the total supernova energy needs to be absorbed by the halo. Is this reasonable? 
Let us assume a kinetic mixing parameter e ~ 10~ 9 , so that around half of type 
II supernova energy is converted into e', e', 7' emitted from the core initially with 
energies ~ MeV. The huge number of energetic e',e', 7' particles injected into the 
region [(~ 1 pc) 3 ] around ordinary supernova will rapidly cool, ultimately converting 
most of their energy into mirror photons. The energy spectrum of these mirror 
photons is naturally very difficult to predict but it could have some resemblance to 
the 7 spectrum of ordinary Gamma Ray Bursts (GRB's). Recall GRB's feature a 
fairly wide spectrum of energies with mean ~ 700 keV with a few percent of energy 
radiated below around 10 keV. In any case, these mirror photons will then heat 
the mirror particle halo, potentially supplying the energy lost from the halo due to 
radiative cooling. 

It has been argued previously [33, 37] that this 7' energy cannot be transferred 
to the halo via elastic (Thomson) scattering off free e' in the plasma. The Thomson 
cross-section is at least an order of magnitude too small. Thus, if the halo contains 
only H' and He' components, then it is hard to see how enough energy can be ab- 
sorbed by the halo to replace the energy lost due to radiative cooling. However if 
the halo contains mirror metal components then the situation is much more promis- 
ing. The heavy metal components are not fully ionized but can have their atomic 
K-shells filled. The photoionization cross-section is many orders of magnitude larger 
than the Thomson cross-section and even a small metal component can make the 
halo optically thick, at least for a range of 7' energies. Once the energetic K-shell e' 
is ejected from the ion, it will interact with the free e' and the ions in the vicinity 
(typically ~ pc) and thermalize. 

The total photoelectric cross-section 2 of a mirror element with atomic number, 
Z, is given by (see e.g. [38]): 



a PE (E Y ) = »^" a 6 Z 5 



16 ^ «-rmel r/2 for£ 7 ,»J (4) 



e 



.Ey 



where I is the e' binding energy and g = 1 or 2 counts the number of K-shell mirror 
electrons present. Evidently, the photoelectric cross-section decreases with mirror 
photon energy like (_Ey)~ 7 / 2 . For Ey near threshold the cross-section has a slightly 
softer behaviour, <jpe oc 1/-Ey and drops abruptly to zero at Ey = I [39]. The 
contribution to the optical depth from such inelastic scattering for 7' propagating 



2 Unless otherwise specified, we use natural units with h = c = 1. 



out from the galactic center is 

tis = J2 2 (Tp E n A 'dr 

~ Y^ 2p r ap E 



m A > 



(5) 



where ^a 1 is the proportion by mass of the mirror metal component, A' (e.g. A' = 
C',0', Si',Fe', ...). The quantities po an d r are the halo central mass density 
and core radius, whose product po r o has been inferred to be roughly constant (i.e. 
independent of galaxy luminosity), with value around 10 2 ' 2 m Q /pc 2 for a Burkert 
profile [40, 41, 42]. 

If we consider just the Fe' component, we find that the optical depth is substan- 
tial, tis ~ 0.1, provided 



E y ~ 30 keV 



Poro 



lO 2 - 2 m /pc 2 



> 



-i 2/7 r t -i 2/7 
KFe> 

0.05 



(6) 



If the ^Fe' component is not too small (~ 0.01) and assuming that the supernova 
7' spectrum peaks above ~ 30 keV and falls sharply at low energy (like a thermal 
spectrum) then the supernova energy being absorbed by the halo will arise (predom- 
inately) from 7' in the energy range 

9 keV ~ E Y ~ 30 keV . (7) 

Even if there are more abundant lighter components, such as O', this may not change 
this picture greatly. Thus, including just the Fe' component, might be sufficient, 
when considering how much supernova energy is being absorbed by the halo. We 
will see in section 5 that the derived galactic scaling properties of spiral galaxies are 
relatively insensitive to the precise details of the supernova 7' spectrum. 

To summarize, ordinary core collapse supernova will produce light mirror par- 
ticles, e', e', 7' from their core with total energy comparable to the neutrino burst 
provided that kinetic mixing of strength e ~ 10 -9 exists. The bulk of this energy is 
expected to be converted into mirror photons, 7', in the region around supernova. 
The details of the resulting 7' energy spectrum are poorly understood, but only the 
part of this spectrum below around 30 keV will be important for heating the halo. 
This heating is achieved by interactions (photoionization) with heavy mirror metal 
components, which occurs because these components retain their K-shell mirror elec- 
trons. Considering just Fe' might be sufficient, as far at the heating of the halo is 
concerned, provided that the proportion of the supernova 7' energy spectrum below 
the Fe' K-shell binding energy, ps 9 keV is small. 

3 Hydrostatic equilibrium 

The halo has two components, a plasma component and a dark disk/compact object 
component. Microlensing observations [43] provide some evidence that the mass of 



the plasma component dominates over that of the compact object component, and 
we henceforth focus on the plasma component. 

The plasma component consists of a set of mirror particles, e' , H' , He', .... What is 
the current chemical composition of the halo? Early Universe cosmology suggests[44] 
that the primordial mirror helium mass fraction is around Y p ps 0.9 for e ~ 10~ 9 , 
with negligible primordial production of mirror metal components. In the first billion 
years or so, substantial mirror star formation and evolution is possible (currently, 
though, the halo is far too hot for much mirror star formation to occur). During this 
early period mirror metals could have been produced reasonably efficiently given that 
mirror stars with large mirror helium mass fraction evolve 10-100 times faster than 
ordinary stars (which have Y p ps 0.25) [45]. Presumably this early epoch of mirror 
star formation is responsible for the current halo metal component (likely at least 
£fv ~ 0.01) required for the halo to absorb enough of the mirror photons produced 
in the hot region around ordinary supernovae 3 . A significant halo metal compo- 
nent is also inferred from the mirror dark matter explanation [26] of the DAM A [27], 
CoGeNT[28] and CRESST-II[29] direct detection experiments. The end result is a 
halo currently composed primarily of He', H', with He' mass fraction around 0.9. 
Additionally there is a small metal fraction (a few percent by mass), which we take 
to be Fe'. An important quantity is the mean mass of the plasma particles: 



m 



^m A ,nA'/^2n A > 



where A' = e', H', He', Fe' and ua 1 is the A' particle number density. For a fully 
ionized plasma, we estimate that m — 1.1 GeV. 

The set of mirror particles, e' , H' , He' , ... in the plasma interact with each other 
via Coulomb scattering. These self-interactions suggest that mirror dark matter 
forms a pressure supported halo. At the present time such a halo would be expected 
to be in hydrostatic equilibrium, where the force of gravity is balanced by the pressure 
gradient. That is, 

— = -p{r)g{r) 



v 2 

u rot 



mn T {r)-^ (9) 



r 



where p(r) = mnT{r), with m being the mean mass of the matter mirror particles 
defined in Eq.(8), nr(r) total dark matter particle number density 4 and g(r) is the 



3 The rampant mirror star formation and evolution in the first billion years or so of galactic evo- 
lution presumably included a large number of mirror supernovae. If the kinetic mixing interaction 
exists, then mirror supernovae should be a source of a large ordinary x-ray photon flux for the same 
reasons that ordinary supernovae are suspected to be a source of a large mirror x-ray photon flux. 
These photons may have been responsible for the reionization of ordinary matter inferred from the 
CMB observations [2]. Importantly, the huge photon flux may not prevent the collapse of ordinary 
matter onto a disk because the ordinary matter with negligible metal component (at that time) 
would absorb a relatively small fraction of the mirror supernovae energy. 

4 If we were to consider also a dark disk/compact object component made of old mirror stars, 
mirror white dwarfs etc, then nr{v) in Eq.(9) would be just the plasma component of the halo. 
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local acceleration due to gravity. Here P(r) = riT(r)T(r) and v ro t(r) is the local 
rotational velocity. 

We shall assume that both the mirror particles and the ordinary baryons are 
distributed with spherically symmetry. Obviously in the central regions of spiral 
galaxies this might not be a good approximation, since the ordinary matter is dis- 
tributed predominately in a disk. However so long as we consider r ~ t\d, then the 
assumption of spherical symmetry could be reasonable. Furthermore, since the typi- 
cal core radius of spirals is inferred to be much larger than the disk scale length, ro, 
it seems probable that the physics responsible for the existence and properties of the 
dark matter core will not depend sensitively on the details of the mass distribution 
at r <C r . 

Anyway, with the assumption of spherical symmetry, g (and hence v^ ot /r) can be 
related to the total mass density, ptotai, via: 

9{r) = ^-f * PtotaidV (10) 



r 2 Jo 



where Gn is Newton's constant. The total mass density, ptotai, can be separated 
into a contribution from ordinary baryons and that due to mirror particles. The 
baryonic contribution of spiral galaxies is approximated by a Freeman disk with 
surface density: 5 



m D _-r/r D Q-p 



2irr 2 D 



e 



where mo is the disk mass and m is the disk scale length. Defining a spherically 
symmetric distribution, pr>, by requiring that the mass within a radius r is the same 
as that of the disk, i.e. J r po^irr^dr' = J r £ 2nr'dr', we have: 

Po(r) = £%-e-""> . (12) 

Studies of spiral galaxies have found that the baryonic mass m D correlates with the 
disk radius, td via[46, 47]: 



lo § I r 5 " I = °- 633 + °- 379 lo £ I in T! D ) + 0.069 
\kpcy \10 iL m Q l 



, / m D 
log 



10 n m m 



(13) 



The distribution of dark matter is something that we would ultimately like to de- 
termine. It might be possible to evolve the halo from an early time and work out 
the dark matter distribution from first principles. In the present work, we will be 
somewhat less ambitious, and assume we know something about the form of the 
dark matter distribution. We will assume that dark matter can be approximated by 
a smooth cored distribution 



p dm = mn T (r) 



Por ° 3 (14) 



(r 2 + ro)(r + r ) 



5 We do not include any other baryonic contribution other than the disk, so that mu represents 
the total baryonic mass of the galaxy. 



where ro, po are the dark matter core radius and central density respectively. Such 
a distribution, known as the Burkert profile, has been suggested by fits to rotation 
curves of spiral galaxies, and other data[48]. 

If Pr>(r) and Pd m {r) are known we can use the hydrostatic equilibrium condition to 
figure out the temperature profile, T(r). To do this, we need a boundary condition. 
Far from the galactic center, i.e. far from heating sources we expect isothermal 
conditions, which motivates dT/dr = at the virial radius r = R V i r . In practice the 
derived temperature profiles have very little dependence on R V i r . To get an idea of 
the typical temperature profiles that we expect for spiral galaxies, we consider three 
examples: (a) a small sized spiral galaxy with disk mass mo = 10 9 m Q and dark 
matter core radius ro = 4 kpc, (b) a medium sized spiral galaxy (~ Milk Way) with 
disk mass mo = 10 n m Q and dark matter core radius r = 12 kpc and (c) a large 
sized spiral galaxy with disk mass mo = 10 12 m Q and dark matter core radius r = 40 
kpc. In each case we take po — lO 2 (m /pc 2 )/r o , consistent with results inferred from 
observations [40, 41, 42]. The results of numerically solving Eq.(9), with the boundary 
condition discussed above are shown in figure 1. The figure indicates that with a 
Burkert dark matter density profile the corresponding temperature profile derived 
from the hydrostatic equilibrium condition smoothly rises towards the central region 
of the galaxy, where it is roughly isothermal. 
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Figure 1: The temperature profile for spiral galaxies for three examples: (a) a small sized 
spiral galaxy with disk mass mo = 10 9 mQ and dark matter core radius r$ = 4 kpc (dashed 
line), (b) a medium sized spiral galaxy (~ Milk Way) with disk mass mn = lO 11 m,0 and 
dark matter core radius vq = 12 kpc (solid line) and (c) a large sized spiral galaxy with 
disk mass mo = 10 12 mQ and dark matter core radius rQ = 40 kpc (dashed-dotted line). 
In each case we take po = lO 2 (m0/pc 2 )/ro. 



4 The ionization state of the halo 



< 



For the temperature range of interest for spiral galaxies, typically 0.01 keV ~ T 
few keV, the plasma is kept ionized by e' collisions. Considering, just the mirror 
helium component for now, the relevant processes are: 



e' + He'° ->• He' + + e' + e' 

J+ , ttV 



e' + He' + -)■ He' z+ + e' + e' (15) 

where He' , He' , He' denote the neutral mirror helium atom, singly charged mir- 
ror helium ion and doubly charged mirror helium ion. The cross-section for these 
processes is known to be reasonably well approximated by the Lotz formula[49] : 



a I = 4.5 x 10~ 14 



ln(E/I) 

EI/eV 2 



cm (16) 



where E > I is the energy of the incident e' and / is the ionization potential. For 
the first process in Eq.(15), / = 24.6 eV while in the second process, / = 54.4 eV. 
We denote the corresponding cross-sections as erf and a\ respectively. 

Opposing ionization are the e' capture processes. The relevant processes for He' 
are: 

e' + He' + ->• He'° + i 
e' + He' 2+ ->■ He' + + 7 ' . (17) 

The cross-section for the capture processes can be approximated by a modified 
Kramers formula [5 0] : 



°c = Y. 



57T a 



5 z. 



eff 



3^3 n 3 E e ,E^ 



where Ey = E e > + e/ ^ 2 — -. For the applications to He', H' and also Fe' ions that we 
will consider, Z e ff = (Zc + Zj)/2, where Zc is the charge of the nuclei and Zj is the 
ionic charge before e' capture [50]. Thus, Z e ff ps 1.5 for the first process in Eq.(17), 
and Z e jf = 2 for the second process. We denote the corresponding cross-sections as 
Uq and a b c respectively. 

The above processes dictate the number density of He' via: 
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n e >n He ,+ (a b jv e >) - n e >n He ,2+ (o b c v el ) (19) 



where the brackets (...) indicate the average over the e' velocity distribution taken 
as a Maxwell-Boltzman distribution: 
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(a b c v e ,) = 2J-— (-) / c'ce-^'l 1 E e ,dE e , . (20) 



In a steady state situation, we have dn He a+ jdt = and thus 

T?He' _ n He' 2 + (gfoe') / , x 

K 2 = = T~b \ • W 

Similarly for the processes affecting dn He ,o/dt, and for the corresponding process for 
mirror hydrogen: 

R He' = n He'+ = K?V) 
1 ~ n He'° (°C V e') 

H , _ n H ,+ (gve,) 

K i = = 7 r K zz ) 

n H ,o \crcVe>) 

where 07 and Oc are the relevant cross-sections for the mirror hydrogen process. The 
H' ionization cross-section is given by Eq.(16) with I = 13.6 eV, while the capture 
cross-section is given by Eq.(18) with Z e ff = 1. 

With these definitions we can determine the number density of each component 
as a function of one of them, which we choose to be n# e / = n He ,o + n He ,+ + n He a+ : 

( R? e 'R^ \ 
UHe ' 2+ ~ {TTrT^TrWW) nHe ' 

( R? e ' \ 

%e ' + _ V 1 + R? e ' + R¥ e 'R¥ e ' ) nue ' 

n He'° = n He' ~ n He ,+ - n He ,2+ 

n H , + = [yTw] fUHe ' 

n H ,o = fn He > - n H ,+ 

n e > = 2n He ,2+ + n He ,+ + n H ,+ 

n T = (1 + f)n H e> + n e i (23) 

where / = nn'/nne' and tit is the total particle number density. The fraction, /, 
can be related to the He' mass fraction: 

*■' = TT774 • < 24) 

Unless otherwise stated, we take / = 0.4 in our numerical work (which, as we 
already discussed at the beginning of the previous section, is suggested from early 
Universe cosmology). The quantities, R^^ and R^ depend only on the temperature. 
It is straightforward to compute the He' ionization fractions, F^ e = n He ,o/n He >, 
Ff Ie = n He ,+ /nHe', Fi 16 — n He ,2 +/ n He' and also the H' ionization fractions, F$ = 
n H io/riH', Ff = n H ,+ /riH>- We show these results in figure 2 for He' and figure 3 
for H '. Figure 2 indicates that He' is nearly fully ionized for T ~ 10 eV. This is 
substantially below the, / = 54.4 eV ionization energy of He' + , which occurs because 
the capture cross-section is several orders of magnitude smaller than the ionization 
cross-section. Qualitatively similar results arise also for H'. 
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Figure 2: The He' ionization fractions, as a function of the local halo temperature, T. 
Shown are F^ e = n He ,o/riHe' (dashed-dotted line), F( Ie = n He ,+ /njje' (dashed line) and 
F^ 6 = n He ,2+/riHe' (solid line). 
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Figure 3: The H' ionization fractions, as a function of the local halo temperature, T. 
Shown are Fq = n H ,o/n H i (dashed-dotted line) and Ff 1 = n H ,+ /tih> (solid line). 
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In addition to the pure H', He' halo, we will consider a small metal component, 
which we take as Fe' with total number density, n Fe i. We denote the number density 
of completely ionized Fe' as n Fe >** and Fe' with 1 K-shell e' as n Fe r* . The ionization 
energy of the bound mirror electron in Fe'* is 9.3 keV and if both K-shell mirror 
electrons are present, the binding energy is 8.8 keV[51]. In figure 4 we show the 



computed ionization fractions, F[ 



Fe' 



n Fe '*/n Fe 



and F^ e 



n Fe >**/n Fe 



versus 



temperature. Figure 4 indicates that Fe' is nearly fully ionized until the temperature 
drops below around 20 keV i.e. somewhat above the ionization energy of the K- 
shell bound mirror electrons. [This occurs because the Fe'* capture cross-section is 
somewhat larger than the ionization cross-section.] For temperatures below around 
2 keV, greater than 99 percent of the Fe' ions typically has both atomic K-shell 
states filled. 
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Figure 4: The Fe' ionization fractions, Ff e = n Fe i* /n Fe i (dashed line), F^ e 
(solid line) as a function of the local halo temperature, T. 
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5 Halo energetics and scaling relations 

Studies of spiral galaxies have found that the baryonic parameters, mD,rD,L r sat- 
isfy two approximate relations. One of these relations was given in the previous 
section, Eq.(13), the other relates the r*-band luminosity, L r to mo (see e.g. [52] 
and references there- in): 

/ r \ L3 
m ° - ' I K * (25) 



10 n m o \Lf w ) 

where Lf w k2x 10 10 L q is the r*-band luminosity of the Milky Way. With these 
relations, the baryonic parameters of spirals are (roughly) specified by a single pa- 
rameter which can be taken as one of mD,fD or L r . If we assume that the dark 
matter is distributed via the Burkert profile, Eq.(14), then we have two further pa- 
rameters, po and tq. Our aim is to derive constraints on these two dark matter 
parameters from dynamical considerations. The derived constraints can then be 
compared with 'empirical' relations derived from observations of rotation curves in 
spiral galaxies[47, 40, 41, 42, 52, 53, 3]: 

log{-^-\ ~ 2.2 ±0.25 



jn Q pc 



,2.65 



3xW ll m® 




i j_ ( t± ^ 

~ ysxio 11 ™.®) 



2.00 



0.66 + 0.58 logl-^p—) (26) 

where m^ is the halo mass. Spiral galaxies typically have halo's in the mass range, 
10 n m Q ~ rrih ~ lO 13 m0 and baryonic mass in the range, 10 9 m Q ~ mo ~ lO 12 m . 

As discussed in refs.[33, 37] and reviewed in section 2, the energy lost in the halo 
due to dissipative processes might be replaced by supernova energy transported to 
the halo via mirror photons 6 . For this idea to work out, the halo must evolve to a 
state such that the energy being absorbed in each volume element must equate to the 
energy being radiated from the same volume element. Thus, we have a dynamical 
condition: 

a E in a E out , . 



dtdV dtdV 

In the following, we derive approximate formula for the left and right-hand sides of 
the above equation. Subsequently we solve the derived equations numerically. We 
will show that the above dynamical condition can be approximately satisfied with the 



6 One can check that the energy transport due to other processes, such as conduction, is negligible 
in comparison to radiation. This is, of course, due to the fact that the mirror photons have a much 
longer scattering length than mirror electrons. 
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Burkert dark matter profile provided that po and ro satisfy certain relations. These 
relations are then compared with the 'empirical' relations obtained from galactic 
rotation curve (and other) data, discussed above. 

The region around a single supernova will be a source of a huge flux of mirror 
photons with total integrated luminosity of around 10 53 erg, provided photon - mirror 
photon kinetic mixing of strength, e ~ 10~ 9 exists. We define the total (time) average 
mirror photon luminosity due to ordinary core collapse supernovae by 

Lsn — RsnJsn{Esn) (28) 

where Rsn is the galactic supernova rate, fsN is the fraction of supernovae energy 
emitted by mirror particles, e',e', 7', and (E SN ) is the average total energy emit- 
ted per supernova. The quantity L' SN obviously depends on the particular galaxy 
concerned. For the Milky Way galaxy, we have: 



iMW _ / n SN \ / JSN \ / \&SN) \ , , w 1A 44 




L'sn « TT^TTT ^ , 'T"'„ 1.4 x 10" erg/ S . (29) 



To proceed further we will need to parameterize the 7' energy spectrum from each 
supernova source. We assume that the peak of this 7' energy spectrum occurs at 
energies somewhat greater than the K-shell e' atomic binding energy of Fe', ~ 9 
keV. In this case only the low energy part of the spectrum can heat the halo. We 
parameterize this energy spectrum via a power law: 

E J2± = (1+sl) (%.)" hNEsN 



dEy \ E c J \ E c 

= k(E y ) Ci . (30) 

This spectrum has been normalized such that 

E c dN 1 

Ey—— dEy = /snEsn ■ (31) 

u.E/y 

We will consider 1 < C\ < 3 in our numerical work (ci = 2 corresponds to a thermal 
spectrum). Although the spectrum would not be expected to be a power law for 
energies sufficiently high, such details will be unimportant since the halo is optically 
thin at energies E' ~ 30 keV. 

Observations indicate [54] that the supernova rate scales with galactic B-band 
luminosity, Lb, via Rsn oc (Lb) ' , with an uncertainty in the exponent around 0.1. 
With the above definitions, we find for J'sn ~ 0.5 (i.e. for e ~ 10~ 9 ), 

k R ~ 1+Cl ( Lb \ r' MW (W\ 

K RsN ~ JEJ+Z [jMW) L sn (32) 

where Lg W k 2x 10 10 L q is the reference B-band luminosity for the Milk Way. 
Although there is substantial uncertainty in L' SN , possibly as large as an order of 
magnitude, the galactic scaling behaviour of k Rsn should be more certain. 
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Galactic Center 

Figure 5: The geometry. Mirror photons travel a distance d from a supernova source at 
point Q, to heat the halo at a point P. 



The energy spectrum of 7' from a single supernova source is given in Eq.(30). The 
supernovae are distributed throughout the disk. We wish to figure out the (time) 
average 7' flux at a point, P, a distance r from the galactic center, from all supernova 
sources. We assume that the supernova distribution traces the baryonic mass density 
[Eq.(12)]. It then follows that the contribution to the 7' energy flux at P originating 
from a volume element, dV = 2irr' 2 dcos8dr', at a point Q is 

d 2 F(r) K(E Y ) Cl R SN e- T p D 



dEydV 



47rd 2 



m D 



(33) 



where k is defined in Eq.(30), d = y/r 2 + r' 2 — 2rr' cos 9 is the distance of the source 
Q to the point P and r is the optical depth along that path. Summing over all 
contributions we find that the total differential energy flux is given by: 



dF(r) k(E y ) Ci R sn {<*> f p D e~r r' 2 

2d 2 



dEy 



m D 



d cos 6dr' . 



The optical depth r is given by 



J2 n i( r i)o-idy . 

i 

The relevant geometry is shown in figure 5, and we have 

r\ ■- 
cosip - 



(34) 



(35) 



y 2 + r' 2 — 2r'y cos ip 



2r'd 



(36) 



The flux of supernova 7' at a particular point, P, will deposit an energy per unit 
volume per unit time of: 



d?E h 



dF 



dEydtdV dEy 



J2 n i( r ) a i ■ 



(37) 



15 



Thus to calculate the heating at a particular point in the halo we need to deter- 
mine the cross-section and number densities of the various components in the halo 
[denumerated by i, in Eq.(35) and Eq.(37)]. The supernova 7' are assumed to have 
relatively high energies, so that the optical depth is dominated by scattering of 7' off 
bound atomic e'. This is possible because heavy elements, such as Si', Fe', are not 
completely ionized but have their atomic inner shells filled. If the flux of supernova 7' 
are falling below 10 keV, then the interactions with K-shell Fe' mirror electrons are 
likely to be the most important (see discussion in section 2). The total photoelectric 
cross-section is given approximately by Eq.(4) and the Fe' number density is given 
by np e ' = nne' (l + 4) (^~) ( x / ; ) ■ [Also needed are the Fe' ionization fractions 
which can be computed as per section 4.] 

Having discussed the heating rate at the point P in the halo, we now turn our 
attention to the cooling rate at the same point. The cooling rate is expected to 
have contributions from three sources: thermal bremsstrahlung, line emission and 
recombination. We first consider the bremsstrahlung component and comment on 
the line emission and recombination contributions in a moment. 

The rate at which bremsstrahlung energy is radiated per unit volume, per unit 
time is [39]: 



d 2 W 16a 3 {2irT\ 
dtdV 3m e \3m e 



\l/2 
) 5Z [ Z j n j n e'9B\ (3* 



where the index j runs over the mirror ions in the plasma (of charge Zj) and gs is 
the frequency average of the velocity averaged Gaunt factor for free-free emission. 
We take cjb — 1-2, which, as reviewed in ref. [39], should be accurate to within about 
20%. 

In principle the energy radiated at a point, P, can have important contributions 
from line emission and recombination in addition to bremsstrahlung. Such processes 
would depend on the detailed chemical composition of the halo. However these pro- 
cesses may not be as important as naive first thoughts suggest. The energy of the 
radiated mirror photons from line emission and recombination is typically close to 
the halo temperature, T. But the halo is generally expected to be optically thick 
to mirror photons of these energies (rough estimates indicate). Thus, the capacity 
of these processes to directly cool the halo is expected to be greatly diminished. 
Bremsstrahlung, on the other hand, generally produces mirror photons of lower en- 
ergy. The energy spectrum of bremsstrahlung is flat for Ey <^C T. and reduces 
towards zero for increasing T like ~ exp(—Ey/T) (see e.g. [39]). We therefore ex- 
pect that the bremsstrahlung process will cool the halo much more efficiently As a 
rough approximation, we set d d fj^ t = e fj^i where ej is an efficiency factor which 
we set to unity in our numerical work. If the bremsstrahlung process is the dominant 
cooling mechanism then d ^°y ex y/T, that is, an increasing function of T. This might 
explain why the system evolves until dE in = dE out . If a region had dE in > dE out then 
this will make T higher in that region which increases also dE out until dEi n = dE out . 
Similarly if dEi n < dE out then this will make T smaller which decreases dE out until 
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dEi n = dE out . It therefore seems plausible that the system will evolve until the 
dynamical condition, Eq.(27) is satisfied everywhere 7 . 

We are now ready to start solving the equations. The philosophy is that if the 
dark matter density can be parameterized by the form given in Eq.(14) then p an d 
tq can be determined by demanding that d dtd y — d d f d y t for each volume element. To 
quantify how well dE in and dE out match, we introduce the quantity, A: 



1 rllr-D 

A ( r O,Po) = -T7— \ 
10r D Jr D 



d Ej n 
dtdV 



d E ou t 
dtdV 



dr . (39) 



The quantity A can be computed numerically via a fortran code. The adopted pro- 
cedure is to input the dark matter density profile and the baryonic density and work 
out the temperature profile using the hydrostatic equilibrium condition as described 
in section 3. Once the temperature profile is known, one can work out the ionization 
state of the halo via the equations described in section 4. Armed with this informa- 
tion, one can proceed to work out Eqs. (34,37) and Eq.(38), the latter we equate to 
d d f d y t ■ The end result is the trivial integration, Eq.(39), to obtain A as a function 
of the input parameters, po, r as well as baryonic parameter, mp [r^ and L are 
obtained from Eq.(13) and Eq.(25)]. 

We start by examining how compatible the Burkert dark matter profile is with the 
d d ^ d y = d d f d °y condition. Define A^ in (r ) as the quantity A(r ,p ) minimized with 
respect to variations in p . Similarly, A^ in (p ) is defined by minimizing A(r ,po) 
with respect to variations in r . In figure 6a we plot A^ in (r ) versus r and in figure 
6b we plot A^ n (p ) versus po. These figures are for an example point with tud = 
10 11 m Q and reference parameters, / = riH> AVhv = 0.4, £p e / = 0.02, L' SN = 2 x 10 45 
erg/s, E c = 50 keV. 



7 Naturally the system is a complicated one and other feedback mechanisms can also be impor- 
tant. For instance, a mismatch of dEi n and dE out can also cause expansion or contraction of the 
halo, which in turn can affect the ordinary star formation rate and thereby readjust dEi n . Such a 
feedback mechanism may also help regulate the star formation rate as suggested by observations 
[see ref.[55] and references there- in for relevant discussions]. 



17 



0.25 




0.15 



0.05 



r (kpc) 



Figure 6a: A^ in (defined in text) versus the core radius, r^ for mj) = 10 11 m©. Plotted 
are various values of c\ which parameterize the hardness of the mirror photon supernova 
spectrum: c\ = 1 (dashed line), ci = 2 (solid line) and c\ = 3 (dashed-dotted line). 
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Figure 6b: Same as for figure 6a, except A^ in is plotted versus po- 
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Figure 6c: dtd ft / dt ^y l versus r/rp for the same example as figure 6a,b. The value of ro, po 
taken are the one's which minimize A. Three curves (almost indistinguishable) correspond 
to the same three c\ values of figure 6a,b. 



The figures show that A has a minimum value of around 0.1. That is 



d Ej, 
' dtdV 



d d ^ d °y to within around 10%. In figure 6c, we plot ^^r/ ^y* (evaluated at the 
Po, ''"o values which minimize A) versus galactic radius r jr^. This is done for the 
same example as figure 6a,b. This figure demonstrates that the Burkert profile 
does a very good job at minimizing A with only small deviations near the galactic 
center, which are unlikely to be important since ordinary baryons dominate the mass 
density there (see e.g. [3] and references there-in). We have found that varying m^ 
and variations in the other parameters give similar results. This demonstrates that 
the Burkert profile is (roughly) compatible with the condition, Eq.(27). Furthermore, 
for each value of m^ we can estimate the values of r and p by minimizing A, as 
we have just done for the particular example with mp = 10 n mQ. This leads to 
relations connecting the baryonic parameters (mo, )"d, L) with the dark matter 
parameters, po,r . Although we expect to derive only two independent relations, 
we explore our results by considering three plots. In figure 7a we plot L versus tq, 
in figure 7b we plot po versus mp and in figure 7c we plot por versus L. In each 
case we consider various values of c\ which, recall, parameterizes the hardness of the 
supernova 7' spectrum. Other parameters fixed are / = nn'/nHe' = 0.4, £p e / = 0.02. 
The quantity k Rsn is given in Eq.(32) with E c = 50 keV and L' SN = 2 x 10 45 erg/s. 
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Figure 7a: Derived galaxy luminosity versus core radius tq for various values of c\ which 
parameterize the hardness of the supernova mirror photon spectrum. Plotted are: c\ = \ 
(dashed line), c\ = 2 (solid line) and c\ = 3 (dashed-dotted line). Also shown (thick solid 
line) is the corresponding empirical galactic scaling relation obtained from Eq.(26). 
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Figure 7b: Derived dark matter central density, po, versus baryonic mass, mp for various 
values of c\. Plotted are: c\ = 1 (dashed line), c± = 2 (solid line) and c\ = 3 (dashed- 
dotted line). Also shown (thick solid line) is the corresponding empirical galactic scaling 
relation obtained from Eq.(25) and Eq.(26). 
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Figure 7c: Derived po?"o as a function of galactic luminosity, Lb for various values of c\. 
Plotted are: c\ = 1 (dashed line), c\ = 2 (solid line) and c\ = 3 (dashed-dotted line). Also 
shown is po r o = 120 777,0/pc 2 , close to the central value of the 'empirical' relation Eq.(26) 
(thick solid line). 



The figures demonstrate that the derived relations for po) r o are compatible with 
the rough 'empirical' scaling relations given in Eq.(26) 8 . The level of agreement 
seems to be very nontrivial. The only parameter adjusted was L' SN which was set 
so that p r had a value ~ 10 2 m Q /pc 2 for m D = lO n m , C\ = 2. 

To check the robustness of these results with respect to reasonable parameter 
variations we have varied £p e /, V SN . [Note that changing E c has the same effect as 
changing L' SN and is therefore not considered.] In figure 8 [figure 9] we show the 
effect of varying L' SN [£i? e '], with the other parameters unchanged. These figures 
demonstrate that the variation of each of these parameters by an order of magnitude 
around our reference values does not greatly modify the L versus r relation (figure 
7a). They also show that po r o is still constant but the value of the constant is 
modified somewhat. We have also found that the above conclusions hold also when 
/ = riH'/nHe', is changed. Thus, the scaling properties demonstrated in figures 7 
remain valid even when parameters are varied. Figures 7,8,9 indicate that the effect 
of varying ci, ^Fe' and V ' SN on p r can be roughly approximated by: 



Poro 



C,Fe' 

0.02 



0.7 r 



V 



MW 

SN 



-,0.7 



2 x 10 45 erg/s 



10 2 m /pc z 



(40) 



8 We have neglected any possible scaling difference between the galactic r* band luminosity, L r 
and B-band luminosity, Lb- 
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On the other hand the L versus tq relation is remarkably insensitive to variations of 
the parameters. Instead of fixing the baryonic parameters via the relations, Eq.(13) 
and Eq.(25), we can consider independent variations ofm^, L and ro- I n doing so, 
we find that r is mainly set by the parameter r^. That is, it is primarily the disk 
scale length that sets the scale for the dark matter core radius. This suggests that the 
primary relation for r is one in terms of ro (not L). Fixing the baryonic parameters 
as per Eq.(13) and Eq.(25), and minimizing A we obtain the approximate numerical 
result: 



r « 3.0 



Td_ 
kpc 



1.1 



kpc . 



(41) 



As with our other results, this is valid over the considered baryonic mass range of 
spirals lO 9 m ~ vtid ~ lO 12 m , and assumed the Burkert profile. 

These results are all very interesting, and among other things, support the premise 
that kinetic mixing is likely close to e = 10~ 9 . Such a value has already been 
identified as a region of interest from the analysis [26] of direct detection experiments, 
such as DAM A [27]. The latter can be explained with e^/^Fv ~ 2 x 10~ 10 . Observe 
that since L' SN oc e 2 [Eq.(3)], Eq.(40) suggests 



Po^o 



^ 



Fe' 



5 X 10" 10 



1.4 



10 2 m /pc 2 



(42) 



This demonstrates the compatibility of the galactic scaling relations with results 
from the direct detection experiments. 
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Figure 8a: Derived galaxy luminosity versus core radius rg for c\ 
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Figure 8b: Derived dark matter density, po, versus baryonic mass, mo for c\ = 2 and 

i r r /MW T iMW 

various values of L SN : L SN 



0.7 x 10 45 erg/s (dashed line) L's„ = 2.0 x 10 45 erg/s 

Uso 
line) is the corresponding empirical galactic scaling relation obtained from Eqs. (25,26). 



(solid line) and L' SN = 6.0 x 10 45 erg/s (dashed-dotted line). Also shown (thick solid 



100 




le+09 



1e+10 

L B [ L sunl 



1e+11 



Figure 8c: Derived po^o as a function of galactic luminosity, Lb for c\ = 2 and various 
values of L'f^. Plotted are L'f^ = 0.7 x 10 45 erg/s (dashed line) L'f^ = 2.0 x 10 45 
erg/s (solid line) and I/gjv = 6-0 x 10 45 erg/s (dashed-dotted line). Also shown is po r o = 
120 itiq/pc 2 , close to the central value of the 'empirical' relation Eq.(26) (thick solid line). 
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Figure 9a: Derived galaxy luminosity versus core radius r^ for c\ = 2 and various values 
of £_F e '. Plotted are ^Fe' = 0.007 (dashed line), t^Fe' = 0.02 (solid line) and £pe' = 0.06 
(dashed-dotted line). Also shown (thick solid line) is the corresponding empirical galactic 
scaling relation obtained from Eq.(26). 
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Figure 9b: Derived dark matter central density, po, versus baryonic mass, mu for c\ = 2 
and various values of £,Fe'- Plotted are £p e / = 0.007 (dashed line), £p e i = 0.02 (solid line) 
and £pe' = 0.06 (dashed-dotted line). Also shown (thick solid line) is the corresponding 
empirical galactic scaling relation obtained from Eqs. (25,26). 
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Figure 9c: Derived po^o as a function of galactic luminosity, Lb for c\ = 2 and various 
values of £i? e '. Plotted are £p e / = 0.007 (dashed line), £p e i = 0.02 (solid line) and £p e / = 
0.06 (dashed-dotted line). Also shown is po r o = 120 m /pc 2 , close to the central value of 
the 'empirical' relation Eq.(26) (thick solid line). 

6 Dwarf spheroidal galaxies 

The results of the previous section are very encouraging. One might be tempted to 
investigate other classes of galaxies, i.e. galaxies beyond spirals, and inquire if they 
are also compatible within this dissipative dark matter picture. A promising class 
of galaxies is the dwarf spheroidal galaxies. These are much smaller than spirals 
and feature luminosity around ~ 10 5 — 10 6 L Q . Observations indicate that these 
galaxies are much more dark matter dominated than spirals. Their low luminosity 
suggests that the halo heating from ordinary supernovae is likely much less than 
typical spirals (~ 4 orders of magnitude less than the Milk Way if the Rsn oc 
(Lb) ' 73 scaling behaviour remains valid). The cooling rate can also be suppressed 
if the halo temperature for dwarf spheroidal galaxies is below around 3 eV. For such 
temperatures figure 2 indicates that most of the mirror helium would be in neutral 
atoms. Only the smaller mirror hydrogen subcomponent is ionized and could thereby 
participate in bremsstrahlung cooling. 
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Figure 10: The temperature profile for a dwarf spheroidal galaxy. Shown is an example 
with dark matter core radius tq = 0.5 kpc and central density po = 10~ 23 g/cm (solid 
line). Also shown is another example with the same dark matter core radius ro = 0.5 kpc 
but slightly higher central density po = 5 x 10~ 24 g/cm (dashed line). 



The halo temperature for dwarf spheroidal galaxies can be computed in the same 
way in which we computed the temperature for spirals, discussed in section 3. Con- 
sider two examples with (a) dark matter core radius r = 0.5 kpc and central density 
Po = 10~ 23 g/cm 3 and (b) r = 0.5 kpc and central density po = 5 x 10 -24 g/cm . 
Both examples are consistent with the observations [56]. The results of numerically 
solving Eq.(9) are shown in figure 10. When computing the temperature profile for 
dwarf spheroidal galaxies, we have neglected the baryon component completely as 
far as its contribution to the local acceleration, g. This is reasonable because the 
baryonic mass is estimated to be only a few percent of the total mass for these 
galaxies[56]. Note that for both examples, figure 2 indicates that only a small pro- 
portion of mirror helium will be fully ionized. [Numerically we find that example (a) 
has mean mass parameter given by fh ~ 1.6 GeV, while example (b) has fh ~ 2.3 
GeV.] 

Without assuming any values for p an d r we can try to derive the value for 
these quantities from dynamical considerations by doing an analysis similar to the 
previous section. Since the properties of dwarf spheroidal galaxies are not as well 
constrained as spirals, we consider a 'typical' dwarf galaxy, with parameters, ro = 0.2 
kpc, Lb = 2 x 10 5 L©[56]. Assuming the Rsn oc (L^) 0,73 scaling behaviour remains 
roughly valid down to the luminosity of dwarf spheroidal galaxies, we can minimize 
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A as per section 5 to obtain estimates of po and tq. We do this for the same reference 
parameters as in section 5, i.e. / = 0.4, £ Fe > = 0.02, L' SN = 2 x 10 45 erg/s, E c = 50 
keV and c\ = 2. Performing this analysis, we find p — 10~ 23 ' 14 g/cm 3 and r — 0.57 
kpc [giving p r o — 61 m Q /pc 2 ]. These values are consistent with estimates from 
analysis of rotation curves of typical dwarf spheroidal galaxies [56]. 

7 Elliptical galaxies and clusters 

Elliptical galaxies are another interesting class of objects to think about. These 
galaxies are thought to arise from mergers of two or more large galaxies. An impor- 
tant characteristic of elliptical galaxies is that they are largely devoid of gas and do 
not exhibit significant star formation at the present epoch and thus the current rate 
of ordinary supernovae is very low. In the absence of supernova heating, the mir- 
ror particle plasma is expected to undergo gravitational collapse onto a disk/bulge. 
For large elliptical galaxies, the time scale for this to occur can be long, and even 
today, many such galaxies might not be fully collapsed. In anycase, in the absence 
of a significant heat source it is reasonable to expect the dark matter within ellipti- 
cal galaxies to be flattened out, to some extent, due to dissipative processes. This 
picture is consistent, at least qualitatively, with the observed ellipticity of elliptical 
galaxies, e.g. [57]. That is, for mirror dark matter, and closely related hidden sec- 
tor models, the non-spherical distribution of the dark matter in elliptical galaxies 
can potentially be explained due to dissipation effects (in the absence of which, the 
observed ellipticity could have been used to place stringent limits on dark matter 
self-interactions, as in e.g. [18]). 

Within larger structures, such as clusters of galaxies, mirror dark matter self- 
interactions might lead to important effects. In particular, observations of the bul- 
let cluster have been used to argue for stringent constraints on dark matter self- 
interactions [58]. Recall, the bullet cluster is an example of a system in which a 
collision between two clusters has taken place. Each cluster has three components, 
the galaxies, hot intergalactic gas and then there is the dark matter. When the main 
cluster and subcluster collide, the hot (ordinary matter) gas associated with the two 
colliding clusters appears to be slowed, but not stopped, by interactions. Both the 
galaxies and dark matter components appear to pass through each other. These 
observations pose a potential puzzle for mirror dark matter. Why doesn't the dark 
matter within each cluster slow down due to interactions? Part of the explanation 
could be that the fraction of dark matter existing as hot gas unbound to individ- 
ual galaxies (i.e. intergalactic gas) might be less than the corresponding fraction 
for ordinary matter [59]. Another part of the explanation could be that mirror dark 
matter has self-interactions which are weaker than those of ordinary matter due to 
environmental conditions. In particular, if the mirror particle plasma has higher 
temperature than the ordinary particle plasma. Recall that the mirror electron scat- 
tering cross-section behaves like da/dQ oc l/t> 4 ~ 1/T 2 . The temperature, T, of 
the mirror particle plasma is uncertain, but the hydrostatic equilibrium condition, 
Eq.(9), suggests that Tocm. The higher helium mass fraction in the mirror sector 
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means that m « 1.1 GeV for mirror dark matter, cf. rh ps 0.55 GeV for ordinary 
matter, and thus the temperature of the mirror plasma in clusters might be higher 
than that of the ordinary plasma by a factor of around two or so. If this is indeed 
part of the explanation of the bullet cluster observations, then studies of colliding 
clusters which feature lower temperatures could find an offset between the galaxy 
component and total mass distribution. Interestingly, there are some tentative hints 
in this direction [60]. 

8 Concluding remarks 

We have examined galactic structure within the context of dissipative dark matter 
candidates, focusing on the mirror dark matter model. At first sight, dissipative 
dark matter might seem unlikely, given the inferred approximate sphericity of dark 
matter halo's in spiral galaxies. However, the idea[33] that there exists a heat source 
to counteract the halo cooling due to dissipative interactions and that this heat 
source could be ordinary supernova seems to be possible. At any rate, it is a specific 
idea, and mirror dark matter offers a specific model in which to study it. 

In this article we have provided a detailed numerical analysis of this whole picture. 
Although a number of assumptions are made, and a few corners may have been cut, 
our results are very encouraging. Our analysis indicates that the inferred dark matter 
scaling properties of spiral galaxies, discovered by Salucci and others, are explicable 
within this dark matter framework. Moreover our results are remarkably insensitive 
to many of the unknown parameters of the theory (such as details of the supernova 
energy spectrum, mirror metal mass fraction, etc). Thus, we have reason to be 
confident that this dissipative dark matter picture, and to some extent the specific 
mirror dark matter implementation of it, may well be on the right track in explaining 
galactic structure. That is, we agree with the title of ref.[3], "Dark matter in galaxies: 
leads to its nature" . 
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